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Time-resolved coincidence imaging of photoelectrons and photoions represents the most complete 
experimental measurement of ultrafast excited state dynamics, a multi-dimensional measurement 
for a multi-dimensional problem. Here we present the experimental data from recent coincidence 
imaging experiments, undertaken with the aim of gaining insight into the complex ultrafast excited- 
state dynamics of 1,3-butadiene initiated by absorption of 200 nm light. We discuss photoion and 
photoelectron mappings of increasing dimensionality, and focus particularly on the time-resolved 
photoelectron angular distributions (TRPADs), expected to be a sensitive probe of the electronic 
evolution of the excited state and to provide significant information beyond the time-resolved pho- 
toelectron spectrum (TRPES). Complex temporal behaviour is observed in the TRPADs, revealing 
their sensitivity to the dynamics while also emphasising the difficulty of interpretation of these com- 
plex observables. From the experimental data some details of the wavepacket dynamics are discerned 
relatively directly, and we make some tentative comparisons with existing ab initio calculations in 
order to gain deeper insight into the experimental measurements; finally, we sketch out some con- 
siderations for taking this comparison further in order to bridge the gap between experiment and 
theory. 



I. INTRODUCTION 

Coincidence imaging techniques, in which the full mo- 
mentum vector of both photoelectron and photoion is 
measured, have been growing in popularity and sophis- 
tication over the last 20 years. The COLTRIMS (Cold 
Target Recoil Ion Momentum Spectroscopy) community, 
in particular, has developed significant expertise in rel- 
atively high-energy and multi-coincidence measurements 
PJE], typically (although not exclusively) utilising syn- 
chrotron light sources in photoionization studies. The 
original aim of COLTRIMS was application to many- 
body collision dynamics, via kinematically complete mea- 
surements of collision systems [4], although the tech- 
nique has since been applied to studies as diverse as pho- 
toelectron diffraction [6 , probing entanglement [7j and 
the investigation of tunnel ionization [HJ |9], as well as 
extensive studies of photoelectron angular distributions 
[TO] [11] . Flat-field and VMI based coincidence imaging 
experiments have also been used by a handful of groups 
[T2HT8] . usually with a focus on lower-energy processes 
as appropriate to experiments based around table-top 
laser sources. The first demonstration of femtosecond 
time-resolved coincidence imaging to study photochem- 
ical processes was over a decade ago [12] [19] , and this 
type of imaging measurement provides the fullest experi- 
mental dataset possible, which can be considered as a 7D 
measurement, or even 8D if one considers the fragment 



* paul.hockett@nrc.ca 



mass spectrum as a distinct observable to the fragment 
velocity distributions; time-resolved coincidence imaging 
therefore provides the best chance of elucidating com- 
plicated, multi-dimensional, excited state dynamics from 
experimental measurements. 

Despite the potential of coincidence imaging, the tech- 
nique has had only a small impact thus far to time- 
resolved measurements generally, and more specifically 
to measurements utilising UV sources [11] . The difficulty 
of applying coincidence imaging in time-resolved, UV 
pump-probe experiments is partly due to the limitation 
of single-particle counting techniques - with consequent 
requirements for long experimental runs and long-term 
experimental stability - which makes time-resolved ex- 
periments particularly challenging; additionally there is 
the inherent difficulty of producing and controlling short- 
pulse UV light. A particular issue is the minimisation of 
background signal from scattered light, which becomes 
a problem on a per photon basis once photon energies 
are above the work function of the materials used in the 
spectrometer [20H23] . The benefit of UV wavelengths is 
that, for many small molecules, the photon energies are 
sufficient for 1-photon pump, 1-photon probe experimen- 
tal schemes, which are ideal for the study of the dynam- 
ics of electronically excited states of molecules. In these 
kind of schemes the laser intensities can be kept low and 
well within the perturbative regime (<C 10 12 Wcm -2 ), 
and the observables take their simplest form. The few 
successful studies to date [12] [16] [I8j |24] have begun to 
explore the power of the technique, but much work re- 
mains to be done. 



2 



The observables provided in a full imaging study pro- 
vide additional information beyond the energy-time map- 
ping of a, now routine, 2D photoelectron or photoion 
measurement, which provide time-resolved photoelectron 
spectra (TRPES) or mass spectra (TRMS) respectively. 
In particular the time-resolved photoelectron angular dis- 
tributions (TRPADs) are sensitive to the electronic struc- 
ture of the ionizing state, so are expected to reveal subtle 
details of the non-adiabatic electronic dynamics [25H29] . 
In the simplest case one might concoct, that of passage 
through a conical intersection (CI) leading to a change in 
the electronic symmetry of the excited state, the PAD is 
expected to change reflecting the non-adiabatic dynam- 
ics and map principally the electronic part of the dy- 
namics [27, 30 j. This expectation can be contrasted with 
the time-resolved photoelectron spectrum, which can be 
considered to be an observable dominated by vibrational 
motions [30], but may additionally map some aspects of 
the non-adiabatic electronic dynamics cleanly depending 
on the ionization correlations [3TH33] . 

However, in the case of large amplitude motions on a 
single electronic state (adiabatic dynamics) changes in 
the PAD are also expected and, furthermore, the PAD is 
energy dependent, so changes to the vertical ionization 
potential (IP) as a function of nuclear coordinates will 
also couple into the form of the observed PAD. These 
factors make the mapping of the dynamics onto the TR- 
PADs non-trivial to understand at both a qualitative 
and quantitative level, and obviate the (relatively) sim- 
ple picture that the TRPADs map only the electronic 
dynamics for all but the simplest of cases; recent compu- 
tational studies of excited state dynamics in NO2 [29] 
have illustrated the response of the TRPES and TR- 
PADs to complicated excited-state dynamics, demon- 
strating the richness of the observable while also sug- 
gesting the difficulty of obtaining detailed insights into 
molecular dynamics via a purely experimental approach, 
with no a priori knowledge of the underlying dynamics. 
Of particular note is the non-isomorphic nature of the nu- 
clear configuration and observable mapping spaces: the 
wavepacket motion on the excited- state, which includes 
dispersion, bifurcation, interferences and other complex 
quantum- mechanical behaviours, does not allow for a 
direct mapping of a given dimension in nuclear coordi- 
nate space onto a given dimension of the observable (e.g. 
photoelectron energy, anisotropy parameter), although 
such mappings may be possible in low dimensionality 
problems such as vibrational wavepackets in diatomics 
[25| [34H36] . However unsurprising this conclusion is, the 
knowledge gap between the observable and the underly- 
ing wavepacket dynamics is often overlooked or ignored 
in treatments of time-resolved measurements. 

Despite these difficulties, the TRPADs provide an ad- 
ditional observable - therefore more information - than 
the TRPES alone, and additional dimensionality to the 
dataset. One demonstration of the utility of this higher 
information content has been the interpretation of exper- 
imental measurements via qualitative/semi-quantitative 



modelling of TRPADs, which can provide insight into the 
mapping of the excited state wavepacket to the observ- 
ables without the need for full ab initio treatments, and 
yield deeper insight into the molecular dynamics than the 
TRPES alone p~6j [37] . Many studies based on velocity 
map imaging (VMI) measurements have also illustrated 
the utility of TRPADs at a phenomenological level (in no 
small part because PADs come "for free" with the tech- 
nique) and a recent review article has surveyed much of 
the work by such "users" of photoelectron angular distri- 
butions [TTj . 

The ultrafast dynamics of 1,3-butadiene (C^Hq) have 
been studied in considerable detail experimentally and 
theoretically, most recently by Boguslavskiy et. al. [38] 
and Levine et. al. [39] (for a more detailed overview of 
the butadiene literature to date see ref. [39]). In the ex- 
perimental photoelectron study, butadiene was excited to 
the bright l 1 !^ state, with UV radiation around 216 nm 
(5.74 eV), and probed via ionization with a time-delayed 
UV pulse at 266 nm (4.66 eV). The computational stud- 
ies, based around an ab initio multiple spawning (AIMS) 
methodology, focussed on describing the excited state dy- 
namics of population on the same bright state and with 
similar energy to the experimental caseQ The conclu- 
sions from these complementary experimental and theo- 
retical studies are in accord, and point to rapid and com- 
plicated dynamics involving fast motion on the initially 
populated bright l 1 ^ state (historically termed the S2 
state, as it lies higher in energy in the Franck-Condon 
region), with twisting about the carbon backbone and 
out-of-plane bending motions leading to highly distorted 
geometries (relative to the planar ground state) on <40 fs 
timescales [39] . At least two minimum energy conical in- 
tersections (CIs), coupling S2 to Si (the optically dark 
2 x A g state) and three CIs coupling Si to So, were found 
to play important roles in the relaxation of the excited 
state. In a wavepacket picture, the initial dynamics from 
the Franck-Condon region to the first CI would corre- 
spond to rapid passage down steep gradients, with lit- 
tle dispersion of the wavepacket along other coordinates. 
Once on Si, the topology would cause the wavepacket 
to split, with parts heading towards each CI, and there 
would be the possibility of more complex dynamical be- 
haviour. Away from the CIs, the non-adiabatic cou- 
pling of the S2 and Si states is strong over large regions 
of the nuclear configuration hyperspace, with significant 
interaction between the states even at the equilibrium 
geometry - for example, the "dark" 2 1 A g state is ac- 
tually found to carry non-negligible oscillator strength 
due to non-adiabatic coupling with the bright state [39] . 
The strongly coupled nature of these states and short 



1 In this case neither the experimental pump bandwidth, nor the 
probe pulse, were included in the simulations. More recently the 
AIMS methodology has been expanded to explicitly include the 
probe process, although this has so far only been applied to the 
photoionization of ethylene with VUV pulses |40| . 
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timescales involved suggests that butadiene can be con- 
sidered as something of a prototypical, perhaps even lim- 
iting or pathological, case for rapid dynamics; as such 
butadiene is a good exemplar of a system where the infor- 
mation available from frequency resolved measurements 
is very limited^ but the speed and complexity of the 
dynamics make time-resolved measurements technically 
demanding and difficult to interpret. 

Experimental measurements have probed the projec- 
tion of these dynamics onto the time-resolved ion yields 
[43| [44] and time-resolved photoelectron spectrum [38] . 
In the former case fast time-constants, <50 fs, were deter- 
mined for the decay of the parent ion signal, and delayed 
onset of fragmentation was also observed [43] ; in the lat- 
ter case, as well as a rapid decay of the photoelectron 
yield, a fast shift in the vertical ionization potential and 
a photoelectron band with little structure were observed. 
Recent work, in which the TRPES was measured in coin- 
cidence with the mass spectrum [38], has demonstrated 
the separation of the photoelectron spectra correlated 
with the S2 and 5i states due to the ionization correla- 
tions with Dq and Di ion states respectively; the Do state 
is stable, yielding parent ion signal only, while the D± 
state fragments [45H47] , thus producing a very different 
mass spectrum. In this way coincidence measurements 
have been able to help disentangle the TRPES data by 
providing complementary, correlated observables. The 
same photofragment coincidence technique has also been 
applied to strong-field ionization of butadiene in order 
to probe multi-electron effects [48]. Measurement of the 
TRPES and time-resolved mass spectrum (TRMS) in co- 
incidence provides a 3D dataset, with the additional di- 
mension of ion time-of-flight and, depending on the de- 
tails of the measurement, may also provide information 
on the kinetic energy release of the fragments. One might 
hope, therefore, that a full, 7D, coincidence imaging ex- 
periment provides a rich enough dataset to discern some 
specific, possibly mechanistic, details of excited-state dy- 
namics in molecular systems, even in the case of very fast 
dynamics exemplified by butadiene. 

In order to explore some of these issues we present 
here our recent experimental results, focussing on a pre- 
sentation of a full dataset to illustrate the richness of 
the 7D coincidence measurements. Along with the data, 
we include a brief description of our coincidence imaging 
apparatus, and a qualitative analysis of the data, with 
a focus on the TRPADs. In future publications, aspects 
of our apparatus will be discussed in fuller detail, and 
a more detailed interpretation of the experimental data 



2 For instance ref. discusses the experimental absorption spec- 
tra, and the lack of information obtainable from the very diffuse 
bands observed (with just three broad features visible over the 
44000 - 51000 cm -1 region studied); ref. [42] provides details 
of ab initio calculations of the optical spectra which require the 
inclusion of a phenomenological dephasing constant to match the 
experimental data, again indicating a gap in the understanding 
of the radiationless relaxation of the excited state. 
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FIG. 1. Spectra and autocorrelation traces for the (a) 266 nm 
probe pulses, (b) 200 nm pump pulses. Spectra show central 
wavelength and FWHM for the Gaussian fit. Autocorrelation 
traces show txc (the FWHM of the Gaussian fit), and the 
corresponding pulse r p duration assuming a Gaussian enve- 
lope. 



and comparison with recent ab initio dynamics calcula- 
tions will be made. 



II. EXPERIMENTAL 

Time-resolved pump-probe measurements were carried 
out with sub-40 fs UV pulses, \ pU mp =200 nm (6.20 eV) 
and Xprobe =266 nm (4.66 eV). Butadiene (1% in He) was 
introduced to the interaction region via a 1 kHz pulsed 
valve (Even-Lavie [49], 150 /xm diameter conical nozzle) 
with a stagnation pressure of ~5 bar. Full 7D coinci- 
dence measurements were performed on a coincidence 
imaging spectrometer (CIS). Details of the apparatus are 
provided in the following sections. 



A. Optical set-up 

Short pulse infrared light (A=800 nm, 35 fs, 1 kHz 
repetition rate) was generated by a standard titanium- 
sapphire based regenerative amplifier system, followed 
by a single pass amplifier stage (Coherent Legend Elite 
Duo). Approximately 700 jiJ of the output was used 
to pump a 3 rd (3a;) and 4 th (4a;) harmonic generation 
scheme, based on sum- frequency mixing in thin BBO 
crystals. 

Calcium fluoride prism pairs were used to compress 
the generated UV pulses, and compensate for the disper- 
sion of transmissive optics in the beam paths (harmonic 
separation and recombination optics, A/2 plate for 3 rd 
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FIG. 2. Schematic of the coincidence imaging spectrometer, 
showing key aspects of the design. The inset shows details 
of the interaction region and electron flight region, including 
electron and ion trajectories; dashed lines represent the initial 
velocity vectors and the solid lines the paths in the guiding 
fields. Dimensions shown are in mm, and the axis definitions 
shown are used throughout this work. 



harmonic, experimental chamber window, propagation in 
air). The output pulses were measured directly with an 
autocorrelator based on 2-photon absorption [50 , and 
the pulse durations in the experimental chamber were 
additionally confirmed via the cross-correlation feature 
obtained from fitting the TRPES data. Output pulses 
on the order of 35 ± 3 fs were measured. Typical au- 
tocorrelation traces, along with pulse spectra measured 
with a UV spectrometer (Ocean Optics Maya Pro) are 
shown in figure [I] 

Control over pump-probe delay was achieved via a 
high-precision linear motor stage (Newport XML210), 
for the data reported herein steps of 10 fs were used. 
The beams were recombined in a collinear geometry, sent 
through a spatial filter to clean the mode and increase the 
beam diameters by a factor of 2, and loosely focussed 
into the interaction region via an /= 1 m focussing mir- 
ror. Beam diameters of the UV at the focussing mirror 
were 4 mm at 3uj and 2.6 mm at 4cu. Focal spot sizes 
in the interaction region were estimated to be ~85 /im 
at 3uj and ~100 /im at Auo. To avoid multi-photon ef- 
fects, and minimise scattered light signal from the Auj 
light, relatively low pulse energies of ~ 200 n J at 3u and 
~ 10 n J at 4co> were used, corresponding to peak inten- 
sities of ~1.9xlO n Wcm~ 2 and -7xl0 9 Wcm" 2 respec- 
tively. 



B. Coincidence imaging spectrometer 

The coincidence imaging spectrometer (CIS) is illus- 
trated in figure [2] An overview is provided in the follow- 
ing, and a more detailed description may be found in ref. 



[20] . The imaging of both electrons and ions is based on 
a flat-field, Wiley- McLaren type geometry. A flat-field 
design was chosen to allow for a large turn-around time 
spread of photoelectrons along the time-of-flight (ToF) 
axis (labelled as the z-axis in figure [5]), as compared with 
a VMI type configuration, and also to provide the most 
direct mapping of (x, ToF) data to initial velocity vec- 
tor, thus allowing for a simple calibration and data back- 
transformation procedure. To create fiat fields, the inter- 
action region is enclosed by grids in the vertical direction; 
the open apertures in the horizontal plane, required to 
admit the supersonic molecular beam and laser beams, 
are removed as far from the interaction centre as possible 
and minimised in spatial extent in order to avoid aber- 
rations and fringing fields. The bottom grid is biased 
(typically +10 to +20 V) relative to the top grid (0 V) in 
order to extract photoelectrons towards the lower flight 
tube. After the photoelectrons have cleared the interac- 
tion region a high voltage pulse (on the order of -0.5 to 
-0.8 kV, <15 ns rise time, supplied by a HVC-1000 pulser 
unit from GPTA) is applied to the top grid to eject the 
photoions to the upper detector. The photoion flight 
tube also contains an acceleration region, and an Einzel 
lens assembly (not used in this work), for further con- 
trol over the ion imaging spectrometer conditions. The 
pulsed valve (not illustrated) is situated in a source cham- 
ber, separated from the CIS chamber by a conical molec- 
ular beam skimmer (Beam Dynamics, 1 mm aperture), 
approximately 45 cm from the interaction region. The 
molecular beam is further skimmed by an adjustable slit, 
discussed below. Partly shown in figure [2] are the baf- 
fle arms, designed to limit the background photoelectron 
signal from scattered light, a particularly severe problem 
from the 200 nm pump pulse. Further details of the scat- 
tered light problem and the baffle system employed here 
will be given in a future publication [23] . 

Both detectors are comprised of a triple stack of 40 mm 
diameter MCPs and delay line anodes (Sensor Sciences 
LLC). Time-of-flight measurement is made via a pick- 
off from the MCP front face, and (x,y) position data 
is obtained from the delay line signals. The ion detec- 
tor is offset slightly (15 mm) along the molecular beam 
axis to compensate for the initial translational velocity 
along this axis. The data read-out and storage is han- 
dled by an electronics rack comprising NIM modules on 
a CAMAC bus backbone, for full specifications see ref. 
[20] , with final output of the processed signals to a PC. 
The timing-resolution, averaged over the full detector 
area, is ~200 ps, limited by the pulse propagation char- 
acteristics over the detector face; the spatial resolution 
is <80 /im for both x and y axes of the electron de- 
tector, and <60 /im for the ion detector [20 . The PC 
also controlled other experimental variables such as the 
pump-probe delay, the shutters installed in the pump and 
probe beam paths, and dwell times at each delay. Data 
was measured with dwell times defined by events (as op- 
posed to laser shots), ensuring good Poissonian statistics 
at all pump-probe time delays even when absolute count 
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rates were low. Pump-only and probe-only signals were 
measured at each delay for dwell times of 1/N less than 
the pump-probe signal, where N was the total number of 
delays set. Each set of N delays defines one experimen- 
tal cycle, and the final dataset presented here comprises 
~ 200 cycles; cycles were kept short (<30 minutes of 
real-time) to minimise the effects of any drifts during the 
course of a cycle. 

Instrument resolution is determined both by the overall 
mechanical design, optimised for low-energy photoelec- 
trons and photoions (<2 eV), the acquisition electronics 
and the fields applied for a given measurement [20 . Ad- 
justment of the extraction fields provides experimental 
control over (x, y) resolution via the choice of the energy 
cut-off (image size), in essence higher extraction fields 
increase the dynamic range of the image at the cost of 
energy resolution, while lower extraction fields maximise 
energy resolution over a reduced dynamic range. This 
consideration is common to VMI and other imaging tech- 
niques [51-54 . In 3D imaging, in common with ID time- 
of-flight measurements, the temporal resolution (z-axis) 
is also affected by the choice of extraction fields due to 
the influence of the fields on the turn-around time-spread 
of the particles, hence temporal spread at the detector 
[55] . Under typical operating conditions, with photoelec- 
trons up to 1 eV extracted with a +13.5 V field, the 
photoelectron energy resolution AE/E was calculated 
from SIMION simulations to be 1% (10 meV) for the 
x-axis, 10% (100 meV) for the y-axis and 3% (30 meV) 
for the ToF axis [20]. Similar figures of 1%, 15% and 
<1% were calculated for the ion (x, y, ToF) resolution at 
0.78 eV with a 200 V extraction field and 600 V accelera- 
tion field [2D] . In practice the ^-resolution was improved 
from these design stage simulations by a reduction of 
the ionization volume along the laser axis, reducing Ay 
from the 2 mm assumed in the resolution figures provided 
above. This reduction was achieved by creating a pseudo- 
ID molecular beam source along the ?/-axis, via the use of 
a piezo- actuated razor blade slit, mounted after the skim- 
mer which separates the source and spectrometer cham- 
bers. In the experiments detailed herein the slit width 
was set at 400 /mi, providing a significant improvement 
on the ?/-axis resolution and allowing for energy slices of 
100 meV to be used throughout the photoelectron data 
analysis over the full 2 eV energy range imaged. 



C. Data analysis &; calibration 

For each event - photoelectron and/or photoion detec- 
tion - a data record is stored, consisting primarily of po- 
sition and time-of-flight (x,?/,ToF), and anode charges 
(QxiQy), f° r the ion and electron. Also stored for each 
data record are various indices allowing correlation of 
the event with pump-probe configuration, pump-probe 
delay, experimental cycle etc. In the current electronics 
configuration only a single electron and/or ion event is 
recorded per laser shot, although in principle multi-hit 



operation of the ion detector is possible and only limited 
by the dead- time of the detector, as defined by the time 
taken for pulses to clear the delay lines (<15 ns). The 
maximum count rate for electrons is therefore the limit- 
ing factor, and is the same as the repetition rate of the 
laser (1 kHz), although for operation in a true coinci- 
dence regime lower count rates may be required to limit 
false coincidence events to a reasonable level [20 , 56] [57]. 
In the experiments presented here the count rates were 
kept to <100 Hz for the dominant parent ion channel, 
which should limit false coincidences to <10% [2D]. It 
is interesting, although unsurprising, to note that higher 
total count rates are permissible in a channel-resolved 
experiment as compared to a single channel case [56] . 
For example, the majority of background electrons from 
scattered light appear at early ToF relative to the main 
signal electrons, so are temporally resolved and do not 
contribute to false coincidences. A more thorough statis- 
tical analysis of the multi-channel case will be presented 
in another publication [57] . 

The data hypercube obtained experimentally may be 
filtered along any or all dimensions in order to examine 
correlations, throw out bad data, retrieve various map- 
pings of the data and so on. For example, charge his- 
tograms allow for the determination of a window of good 
events, defined as single hit events, and the rejection of 
bad events where two or more hits were registered within 
the delay line pulse transit period, such bad events have 
higher (Q x ,Qy) than single hit events. Correlation of 
photoelectrons with photoions of high translational ve- 
locity elimi nates s ignal from ionization of background gas 
(see section III A ) ; correlation with a given mass provides 
fragment-resolved data, and so forth. 

Calibration and backtransformation of the data is per- 
formed via a three- step process: (1) the image centre 
(#o,2/o>2o) * s defined, either by manual inspection, or 
taken as the peak in the (x, y, ToF) histograms; (2) the 
(x, y) position bins are converted to mm positions from 
the image centre, based on a (static) look-up table cal- 
ibration, and converted to velocities (V x ,V y ) by making 
use of the measured time-of-flight; (3) a ToF look-up 
table, based on numerical trajectory calculations which 
account for the voltages applied, is computed and applied 
to convert the ToF data to V z . Because of the flat-field 
arrangement the trajectory calculations are straightfor- 
ward and not computationally demanding. 

Once converted to velocity space, energy and angular 
data can readily be extracted from the dataset. Obtain- 
ing these integrated observables is again just a case of 
creating histograms of the various quantities of interest, 
combined with filtering as described above, to provide 
maps of observables for given regions of the data hyper- 
cube. The histogram of events with energy {E = ^mV 2 ) 
versus time delay (t) provides the 2D map (£", t) - the 
TRPES - with bins of width AE and At. Conversion 
of the data to spherical polar coordinates (E, 0, <p) allows 
binning into volume elements AE sin 0A0A<j) for each de- 
lay t, to create 3D maps of the energy and angular dis- 
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tributions of events. Although the full, quasi-continuous, 
3D distribution may be visualised as a set of isosurfaces 
or projection planes from these maps, the data is perhaps 
represented most tractably as a series of TRPADs at se- 
lected energy and time slices, with angular distribution 
1(0, (j); E, t). For cylindrically symmetric distributions, 
further integration over <p can also be performed leading 
to a reduced form 1(6; E, t). Low event number datasets 
also benefit from this integration because without inte- 
gration or smoothing the 3D maps may be too sparse or 
noisy for further analysis. 

The extracted angular distributions can also be de- 
scribed phenomenologically by @lm parameters, 

1(0, 0; E,t) = J2 PlmY lm (0, 0) (1) 
LM 

where IlmI^^) are spherical harmonics of rank L and 
order M. This distribution is general to any scatter- 
ing system [58] but the expansion is constrained, by the 
symmetry of the experiment, to only certain values of 
L, M [59]; in the case of the pump-probe experiment 
considered here, with total absorption of two photons, 
linearly polarised light and pump and probe polarisation 
vectors parallel, the final angular distributions contain 
terms L = 0, 2, 4 and are cylindrically symmetric (M=0) 
with symmetry axis defined by the laser polarisation. For 
the photoelectrons, the (3lm contain details of the ion- 
ization dynamics, and are complex functions of molecular 
geometry and photoelectron energy [26 , 28, 59 . For the 
photoions, the Plm contain details of the fragmentation 
dynamics (although the expansion is not usually written 
in this exact form, see for example refs. [6QH64] ). 
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1+1 'cut-off (~1.8eV) 

FIG. 3. Electron imaging summary. One quadrant 
of the raw, time-integrated 3D electron data is shown 
as nested isosurfaces. The 2D image planes show the 
(x,y), (y,z) and (x,z) projections of the raw data. For 
the position data, 1 bin w 5 /xm. This is an interac- 
tive figure in some versions of this manuscript; the in- 
teractive version, and source data, is also available at 
|http://dx.doi.org/10.6084/m9.figshare. 106343| 



III. RESULTS & DISCUSSION 

In this section we present a complete dataset of pho- 
toelectron and photoion measurements. An overview of 
the data is first given, providing time-integrated 3D visu- 
alizations of the full datasets. Low dimensionality (TR- 
PES, TRMS) and high dimensionality (TRPADs, corre- 
lated data) maps extracted from the data are then pre- 
sented, followed by a brief discussion of the results in the 
context of gaining insight into the underlying molecular 
dynamics. 



A. Overview - visualizing multi-dimensional 
measurements 

The most direct way to begin considering the data is 
via 3D maps of the electron and ion signals. Here we 
present only the time-integrated data in this form, al- 
though visualization of time-sliced data as 3D maps is 
also a useful technique for qualitative analysis. The 3D 
maps provide direct information on the performance of 



the instrument, a quick check for artefacts, some infor- 
mation on the shape of the data in terms of complexity, 
ion fragment yields, energy spectra and so on; they rep- 
resent a framework within which to begin a more detailed 
analysis of the dataj^] 

The overall 3D photoelectron distribution obtained 
(integrated over t) is shown in figure [3j A few of the 
major features in the dataset are labelled. Of particular 
note is the signal along the laser axis leading to blur- 
ring in this direction. Although the interaction region 
along the laser axis is minimized by use of a pseudo ID 



molecular beam, as discussed in section II B[ ionization of 



background gas along the laser axis cannot be prevented. 
However, in a coincidence measurement this background 
signal can be removed at the analysis stage by filtering 
the electron data for coincidences with parent ions in the 
central ion spot (see figure [4] for the corresponding ion 
distribution). A slight up-down asymmetry is visible in 



Figures [3] are interactive in some versions of this manuscript. 
Interactive versions, and source data, are also available at 
|http://dx doi.org/10. 6084 /m9 . figshare. 1063431 
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FIG. 4. Ion imaging summary (time-integrated), (a) Raw 
3D data (x,y,ToF) for the parent ion region. Raw 3D data, 
sliced along the centre of the distribution, is shown by the 
nested isosurfaces, and 2D projections onto the (x, y), (y,z) 
and (x,z) image planes are also shown. For the position 
data, 1 bin ^5 /xm. Colour mapping shows logio (counts), 
with a dynamic range of 5 orders of magnitude, (b) Time- 
integrated mass spectrum. Obtained by integrating the data 
over (x, y) (background gas signal excluded). The inset shows 
a detailed view of the parent ion region, 50 - 58 a.m.u. (c) 
As (a) but for fragment region centred at 39 a.m.u. This 
is an interactive figure in some versions of this manuscript; 
the interactive version, and source data, is also available at 
|http://dx.doi.org/10.6084/m9.figshare.l06343| 



the images; the photoelectrons in the upper hemisphere 
are those which initially have trajectories away from the 
detector, and therefore spend longer in the interaction 
region as they must be turned-around by the extraction 
fields. The asymmetry is therefore ascribed to a combi- 
nation of possible effects: slight inhomogeneities in the 
extraction fields, especially near the grids, and perturba- 
tion by weak external fields. Such effects could perturb 
the photoelectrons, and would become more significant 
for longer interaction time-scales. A central artefact is 
present around (V x , V y ) = and to long ToFs, the source 
of this artefact is unknown, but such structure could arise 
from the creation of metastable autoioinzing states with 
ns lifetimes (e.g. high- lying Rydberg states [65] [66]). The 
strong probe-only signal, which appears in a narrow en- 
ergy band E < 0.3 eV, is clearly visible near the centre of 
the distribution and peaked along the laser polarisation 
axis. For E > 0.3 eV all electrons are from the pump- 
probe signal and, due to the rapid dynamics resulting in 
a rapid shift of the IP (see figure [6| , appear quite dif- 
fuse in the 3D representation with no obvious angular 
dependence in the t integrated data. 

Figure [4] gives an overview of the complete ion data 
obtained (again integrated over all delays t), showing the 
full 3D map and projections onto 2D planes for (a) the 
parent ion mass region around 54 a.m.u. and (c) the 
fragment region around 39 a.m.u. The central panel (b) 
shows the calibrated mass spectrum obtained by inte- 
grating over (x,y), excluding the background gas signal 
(x < 4000). The 3D ion map provides much more direct 
information than the equivalent mapping of the electron 
data, and from the figure various key features are im- 
mediately visible: (1) the main parent ion feature, lo- 
calised in (x,y,ToF) and intense; (2) satellite features 
to the parent ion, assigned as isotopes to higher mass 
(parent +1 and parent +2) and hydrogen loss channels to 
lower mass (parent- 1, parent-2 and parent-3); (3) weak 
fragment channels appearing at shorter ToF, assigned as 
fragmentation of the parent ion; (4) a stripe of signal to 
the edge of the detector, arising from ionization of back- 
ground gas, i.e. molecules not entrained in the molecular 
beam which have zero net translational velocity, so ap- 
pear along the laser propagation axis. The benefits of 
(x, y) sensitivity to the mass spectrum are immediately 
obvious: the background signal is easily gated out of the 
analysis. Mass resolution of < 1 a.m.u. is readily ob- 
tained under these experimental conditions (ion extrac- 
tion pulse of 550 V, acceleration voltage of 150 V, mass 
range - 230 a.m.u.), and butadiene isotope features are 
cleanly resolved in the mass spectrum. 

Figure 4fa), and the inset to|4jb), shows an expanded 
view of the parent ion and nearby features. Underlying 
the main feature, and centred on it, is a diffuse signal 
with the same ToF. This feature is approximately 3 or- 
ders of magnitude weaker than the main feature. This 
diffuse feature is assigned to a combination of (1) parent 
ions which fragment after extraction from the ionization 
region (i.e. within the flight tube) and (2) parent ion 
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formed from dissociation of clusters within the ionization 
region. In both cases there would be a kinetic energy re- 
lease relative to the direct parent ion signal, which would 
lead to the broadened, isotropic feature observed. The 
dimer peak at 108 a.m.u. was also observed (not shown 
in figure [4^b), normalized intensity 2xl0 -4 ), and dimer-1 
and dimer +1 features were also just visible (normalized 
intensities <5xl0 -5 ), confirming the presence of dimers 
in the molecular beam. In the work discussed herein no 
attempt was made to investigate the diffuse feature fur- 
ther by, for instance, varying the gas mixture or pulsed 
valve timing. The diffuse nature of this feature, relative 
to the main parent-ion feature, means that most of these 
ions can be gated out of the analysis, while those that 
remain under the main feature hardly contribute to the 
total signal. 

Further weak, somewhat diffuse and isotropic features 
are observed at earlier ToF than the parent ion. These 
features are assigned to hydrogen loss channels, resulting 
in the species , and . These channels 

are 2-4 orders of magnitude weaker than the main parent 
ion signal. They show some kinetic energy release, result- 
ing in their diffuse appearance in the imaging data and 
a broadening of the peaks in the (x, y) integrated data, 
and this broadening in the ToF clearly increases with 
the number of H- atoms lost (see inset of figure [4jb)) . 
The isotropic nature of these distributions indicates the 
slow release of the fragment, relative to the timescale of 
molecular rotations (~10s of picoseconds), as would be 
anticipated from the complex dissociation dynamics. 

Figure [4jc) shows the imaging data for the fragment re- 
gion around 39 a.m.u. The main feature is assigned to the 
methyl loss channel, resulting in a C 3 Hf fragment. The 
higher mass fragment is assigned to the same channel for 
the +1 isotope, based on the intensity ratio, but could 
also contain contributions from the lower mass 

channel is assigned to C^H^ • As with the hydrogen loss 
features, these channels show isotropic angular distribu- 
tions, but with a larger spread of kinetic energy release. 
At longer ToFs a diffuse tail is observed. The length of 
the tail - stretching to the parent ion feature - indicates 
a relatively long lifetime for the fragmenting complex, 
because fragmentation must occur in all regions of the 
spectrometer to result in the large smearing out of the 
fragment mass spectrum observed. Therefore the upper 
limit for the timescale of fragmentation is the time taken 
for ion extraction from the interaction and acceleration 
region of the spectrometer, hence ns to /is timescales. 
Conversely, ions which appear at the fragment ToF must 
fragment before extraction, so cannot take longer than a 
few ns to fragment. The large range in timescales here 
suggests that multiple fragmentation pathways may play 
a role. 

Also shown in the mass spectrum, figure |4jb) , is an- 
other fragment region centred at 28 a.m.u. This is as- 
signed as ethylene cation, C2H4 . The satellite peaks, 
following the same pattern as the methyl loss region, are 
assigned to the +1 isotope at 29 a.m.u., and further hy- 
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FIG. 5. Time-resolved ion yields for a selection of the mass 
channels observed, (a) Major ion channels. The parent ion 
signal is scaled down by a factor of 20 for plotting purposes, 
(b) Minor channels. 



drogen loss resulting in C2H3 and C2H2 • This region 
has no diffuse tail, indicating more rapid fragmentation 
than for the methyl loss channels. No lower masses were 
observed in the mass spectrum (down to the mass cut- 
off at 3 a.m.u.), indicating that all lower mass fragments 
were produced as neutrals. 

At this level of representation, a few key aspects of the 
dynamics can be discerned: the electron data appears 
relatively structureless, hence spectrally broad and pos- 
sibly varying rapidly as a function of t; the ion data shows 
several fragmentation products, and some limits to the 
timescales of fragmentation for different channels can be 
intuited; the fragment angular distributions appear near 
isotropic, although quantitative analysis is required to 
make this conclusion definitive. 



B. Low-dimensionality mappings: time-resolved 
mass spectrum and photoelectron spectrum 

The next stage in analysis complexity is low- 
dimensionality mappings: here we consider 2D mappings 
of the dynamics. First, the time-resolved mass spectrum, 
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the mapping of product yields vs. time for each mass 
channel; second, the time-resolved photoelectron spec- 
trum, the mapping of electron yield vs. time for each 
kinetic energy channel. 

Figure [5] shows the time-dependence of a selection of 
the mass channels. To obtain this representation, the 
data was integrated over the (x, y, ToF) coordinates for 
each of the features of interest at each pump-probe delay 
t, and was converted from raw counts to count rates to 
allow for the correct weighting of the pump only and 
probe only signals, which could then be subtracted from 
the full pump-probe signal. 

The parent ion signal is observed to rise with the cross- 
correlation, plateau and then fall with a Gaussian tail. 
The diffuse part of the parent ion signal shows a very 
similar response, as do the isotope peaks. The major hy- 
drogen loss channel, C^H^ , rises at later t, and peaks 
around 20 fs after the parent ion. The minor hydrogen 
loss channels, the parent-2 and parent-3 features shown 
in figure |5|b), both peak at around 30 fs after the par- 
ent ion. Therefore, the data indicates that evolution on 
the excited state of around 20 - 30 fs is required before 
these fragmentation channels are open. Similarly, the 
fragment channels assigned to C^H^ and C2H4 show 
peaks around 20 fs after the parent ion. In all cases the 
peak shapes are similar, with non-Gaussian tails. Assum- 
ing that fragmentation only occurs on the Di surface, the 
timescales here indicate that direct ionization to D\ is 
possible very rapidly, and indicates a significant lowering 
of the vertical IP from the Franck-Condon region. At 
the ground state equilibrium geometry Dq and Di lie at 
9.07 eV and 11.39 eV respectively [6Z|, so are separated 
by ~2.3 eV, with Di lying 0.53 eV above the available 
1+1' photon energy of 10.86 eV. Hence the observed dy- 
namics suggest that the vertical IP to D\ falls by ~0.5 eV 
in ~20 fs, and this inferred drop is very similar to the 
shift observed in the photoelectron signal, discussed be- 
low. Another possibility is that the cross-section for ab- 
sorption of a second probe photon increases dramatically 
over the first 20 fs of the dynamics. Such 1+2' processes 
would provide 15.52 eV of energy, allowing population of 
several higher- lying cationic states [46], either via direct 
2-photon absorption, or sequential absorption via Do. 

The low laser fluences used experimentally suggest that 
all observed signals arise from 1+1' pump-probe pro- 
cesses, conversely the known appearance energies of the 
fragments (>11.3 eV [46j [47]) suggest 1+2' processes, 
may be energetically required for fragmentation to oc- 
cur In order to check more carefully for 1+2' processes 



4 Experimentally these studies could only probe the ground state 
minimum geometry, so do not rule out the appearance of frag- 
ments at lower energies as a function of nuclear coordinates, 
which would only be limited by the asymptotic energy of the 
fragments. However, the energetics of the dissociation pathways, 
including transition states and final products, were also consid- 
ered computationally in ref. (4^, and the theoretical results also 
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FIG. 6. TRPES mapping, (a) Raw counts per (E, t) bin. The 
dashed line follows the tilt of the band maximum, highlighting 
the chirp of the signal, with gradient ^20 fs/eV. (b) As (a), 
but normalized to the maximum counts for each E bin in order 
to show the temporal behaviour at each energy independent 
of total counts. The second dashed line follows the low edge 
of the signal, showing the slight broadening as a function of 
energy, with gradient ^30 fs/eV. 

the time-resolved ion yields were also extracted in coin- 
cidence with photoelectrons of energy 0-1.8 eV. This 
produced essentially the same traces as shown in figure 
|5j except slightly noisier. This does not completely rule 
out 1+2' processes, which could also produce photoelec- 
trons in this energy region by population of cation states 
with internal energies of >13.72 eV, or via the sequen- 
tial process of ionization and subsequent excitation of the 
cation. A careful probe power study would be required 
to firmly answer this question experimentally, but this 
remains for future work. 

In summary, the fragment yields peak around 20 - 30 fs 
after the parent ion, indicating rapid dynamics lead to 
the opening of the observed fragmentation channels. Ad- 
ditional filtering of the data for coincidences with elec- 
trons in the 1-photon ionization region suggested, but 
cannot rigorously confirm, that there is no significant 2- 
photon ionization contribution. In a future publication, 
an extended analysis of the time-resolved ion data will 
be made in order to examine the kinetic energy release 
spectra of each fragment. Combined with the available 
fragment energetics data this information may be suffi- 
cient to accurately determine mechanistic details of the 
various dissociation pathways. 

We next consider the {E, t) mapping of the electron 
data, the TRPES, shown in figure [B| The TRPES shows 
an energetically broad photoelectron band up to the 1- 
photon cut-off at ~1.8 eV. The low-energy cut-off at 
~0.3 eV corresponds to the region of probe-only signal. 
The main feature at ~1.2 eV follows the cross-correlation 



indicate that the observed fragments are not energetically acces- 
sible via 1+1' processes. 
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of the pump and probe pulses (r xc ~60 fs FWHM); 
outside of the cross-correlation region the signal rapidly 
moves to lower energies, with a broad stripe visible in 
the data. This stripe is slightly sloped or chirped (this is 
especially clear in the energy normalized representation 
in figure |6jb)), with the onset time varying as a func- 
tion of energy. The chirp appears close to linear over the 
observed energy region, with a gradient of 20±10 fs/eV. 
The temporal width of the band increases slightly as a 
function of energy (shown by the second dashed line in 
figure [(^b)), and the trailing edge of the band follows a 
gradient of approximately 30±10 fs/eV, although there 
is increased, non-linear, broadening below ^0.4 eV. For 
the linear region the broadening of the band is therefore 
~10 fs/eV; a higher temporal resolution measurement - 
smaller At - would provide a more accurate figure. Apart 
from the slope of the band there is little structure ob- 
served, and most of the population rapidly leaves the 
observation window of the measurement (on the order of 
r^c, ~60 fs). This data is very similar to previous TRPES 
studies [38] , but with shorter r xc and obtained for only a 
small set of delays in order to ensure good statistics for 
the TRPADs. Although the clipping of the data along 
the temporal coordinate makes rigorous determination of 
the lineshapes impossible for the peak of the signal in this 
dataset, based on the previous TRPES data the tempo- 
ral profile of the signal is assumed to be approximately 
Gaussian for most, if not all, energies (see also the line- 
outs in figure [8]), with a small (few percent of the total 
counts) non-Gaussian tail outside of r xc . 

In terms of the dynamics, the delayed onset as a func- 
tion of energy is phenomenologically consistent with the 
picture of fast wavepacket propagation on a steep po- 
tential energy surface. Such motion would map primar- 
ily to the vertical IP (assuming that the excited state 
and ionic ground state potential energy surfaces are not 
topologically identical), hence the kinetic energy of the 
observed photoelectrons. The TRPES provides informa- 
tion on the speed of the IP change, and is consistent with 
the appearance of fragmentation channels with a delayed 
onset which require a significant drop in the vertical IP 
to Di relative to the Franck-Condon region. Conversely, 
wavepacket dynamics which cause little change to the IP 
would be responsible for signal observed in a given en- 
ergy region at long delays, arising from population which 
remains within the observation window of the measure- 
ment (although such population may still trace complex 
trajectories in energy space). Since the signal in all re- 
gions outside of the cross-correlation time-scale is small, 
it is clear that wavepacket motion along these coordinates 
is a very minor contribution to the dynamics. The data 
also shows that the observed signal stays near Gaussian 
for all energy slices, again indicating that there is little 
dispersion of the wavepacket. 

In broad terms, the ion and electron 2D data give some 
insight into the wavepacket dynamics, with the general 
characteristics apparent and some clocking of this mo- 
tion possible. The observed energetic shift in the pho- 



toelectron signal is consistent with the vertical IP drop 
required to access Di, hence observe fragmentation, as 
already inferred from the ion data. In order to assemble 
a more detailed picture we next consider the additional 
information available from the TRPADs. 



C. High-dimensionality mappings: time-resolved 
photoelectron angular distributions 

In order to extract TRPADs the data is calibrated and 



rebinned in polar coordinates as detailed in section |II C 
to give the intensity (counts) per 4D volume element 
A(/>sin<9A<9A£;At, denoted 1(0,0; E, t). The choice of 
binning is, naturally, limited by the experimental time- 
steps and instrument resolution; however coarser binning 
can be used in order to improve the statistics at the loss 
of resolution. For the data presented here At = 10 fs 
and AE = 0.1 eV, and in all cases shown here the data 
was also integrated over <\> which, due to the cylindrical 
symmetry of this experiment, results in no loss of infor- 
mation but does improve the statistics per 9 bin. To ob- 
tain the cleanest possible TRPADs the upper hemisphere 
electrons were discarded at the cost of a factor of two in 
counts. The data was also filtered for coincidences with 
the main parent ion feature. Because the parent ion dom- 
inates the signal, this filtering was not required to obtain 
TRPADs correlated with a given ion channel (although, 
more generally, could be used in this way), but did serve 
to remove all signal from background gas and scattered 
light, resulting in cleaner TRPADs albeit at the cost of 
total counts. Figure [7] shows an example of the TRPADs 
obtained in this way for four energy slices. In this rep- 
resentation, the areas of the TRPADs are normalised to 
unity in order to allow comparison of the form of the 
PADs independent of total counts. The error bars and 
scatter of the data points should therefore be used as a 
guide to the statistical significance of the extracted PADs 
over the various energy and time slices shown. The solid 
lines show a fit to an expansion in spherical harmonics, 
as defined by eqn. [I] 

From the data it is immediately clear that the TR- 
PADs exhibit complex behaviour. The TRPADs change 
rapidly as a function of t, with changes on the order of 
10 - 20 fs apparent. Because the laser pulses used in this 
work were around 35 - 40 fs in duration, there is already 
significant temporal blurring in the measured data. De- 
spite this, the changes are still clear and unambiguous 
in the data. For comparison, TRPADs extracted for the 
probe only background (figure (7ja)) show no significant 
changes temporally beyond the signal noise, and conse- 
quent statistical variation of the fit, signifying that no 
temporal artefacts are present in the raw data or intro- 
duced via the data processing. The scatter in the data is 
worst at the poles of the distribution due to the sin nor- 
malisation factor which serves to amplify noise near the 
poles; conversely the equatorial region shows very little 
variability. There is some asymmetry present in the data 
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FIG. 7. TRPADs, I(6;t,E), extracted from the dataset as detailed in the main text, (a) Probe only data, 0.1 - 0.2 eV. (b) - 
(d) Pump-probe data for 0.7 -0.8 eV, 0.6 - 0.7 eV and 0.5 - 0.6 eV respectively. Data points are shown with statistical error 
bars, solid line shows fits to eqn. [TJwith L = 0, 2, 4 and M = 0. 



which is ascribed to a combination of noise (scatter) and 
detector artefacts/inhomogeneities. For this latter rea- 
son the form of the extracted PADs (as defined by the 
fitted Plm parameters) may not be highly accurate, but 
the results do have high precision and reproducibility as 
shown by figure [T^a), so any relative temporal changes 
observed should be reliable and robust within the statis- 
tical uncertainty. 



One challenge of high-dimensionality datasets is the 
presentation and/or reduction of the data to a more 
tractable form to allow for pattern recognition at a phe- 
nomenological or quantitative level. For the TRPADs 
one can reduce the full dataset to the fitted TRPADs, 
which can then be represented as (E,I{6)) or (t,I(0)) 
surfaces, or maps of Plm{E^€). Figure [8] shows examples 
of the TRPADs represented as (£, 1(0)) surfaces (in polar 
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FIG. 8. Fitted TRPADs, 1(0; t,E), as per figure [7] represented as polar surface plots. White lines show the fits at 10 fs 
intervals, the surface and colour map interpolate between these discrete measurements. Lower panels show the fitted /?2o(£) 
and ^40 (£), as well as the counts, for each energy region. 



space) for the same three energy slices as figure [7^b) -(d), 
and the associated Plm (t) parameters are plotted in the 
insets with error bars. Such representation, as a func- 
tion of energy or time, readily allows for comparison of 
the evolution of the shape of the PADs, although visual 
information about the quality of the raw data and fit fi- 
delity are lost, so care must be exercised when drawing 
conclusions from such maps. For the TRPADs shown 
in figure |8j the comments pertaining to figure [7] can be 
reiterated, namely that the TRPADs for the different en- 
ergy regions are significantly different, and show complex 
temporal evolution. More specifically, for the 0.5 - 0.6 eV 
region, the observed TRPADs show a 4-lobed structure 
(significant /?4,o) at t = —10 fs, which evolves to a 4-lobed 
structure with a different orientation at t = 10 fs, returns 
to near its initial shape at intermediate £, and finally 
shows fast oscillations at t > 60 fs, although these later 
oscillations may not be reliable due to the low statistics 
in this region. The colour mapping on the surface plots 
emphasises the evolution of the signal intensity along the 
equator and at the poles of the TRPADs; the correspond- 
ing (3lm plots show that both /?2,o an d /?4,o change signif- 
icantly as a function of t, with /3 2j o displaying two peaks 
in the range 20 < t < 70 fs. The /^o trace shows a 
minimum at t = 40 fs and a maximum at t — 60 fs, 
which coincide with a mimimum and maximum in the 
/?2,o trace, but at earlier times (t < 40 fs) and later times 
(t > 60 fs) the behaviour does not appear to be directly 



correlated to the /?2,o trace. 

For the 0.6 - 0.7 eV energy slice the picture is quite 
different. For the first time step the PAD again shows a 
4-lobed structure, but with intensity peaked at the poles 
and equator, as opposed to at 45° as per the 0.5 - 0.6 eV 
window. This corresponds to a large and positive /J^o, 
compared with a large and negative (3^o for the lower 
energy slice. The (3^o value goes negative, with small 
magnitude, for the following time slices, and appears to 
show a slight oscillation with minima at t = 0, 40 fs and a 
peak at t = 20 fs; the /?2,o parameter shows a correlated 
oscillation but centred around a mean value of ~ 0.5. 
In the surface plots, this oscillation appears as a slight 
breathing of the TRPADs, with the largest changes at the 
poles. At later times, t > 50 fs, more complex behaviour 
is observed, with a significant beating around the equator 
of the distributions as well as at the poles. 

For the 0.7 - 0.8 eV energy slice the behaviour is again 
different, with much less variability in the observed TR- 
PADs over the peak of the signal. The /?2,o value de- 
creases gradually from the local maxima at t = fs until 
t = 50 fs, then increases gradually to t = 70 fs. The /^o 
trace shows correlated local maxima at t = 0, 70 fs, but 
much more variability between these peaks. As was the 
case for the lower energy slices, the data at long delays 
shows significant scatter and has low statistics, so should 
be treated with care and carefully compared with the 
plots showing angular data points (figure [7]) before draw- 
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ing firm conclusions as to the veracity of the temporal 
evolution observed in this region. However, the fact that 
the observed /?2,o dips at the same time (t = 80 fs) over all 
three energy slices, which were analysed independently, 
suggests this behaviour is genuine and not the result of 
random noise. In terms of the form of the TRPADs, 
the oscillations in the /3l,m describe significant oscilla- 
tions which include large changes to the photoelectron 
flux around the equator of the distributions. Because 
the equatorial region is statistically more reliable than 
the noisier poles this again suggests that these observa- 
tions are genuine. 

Figure [9] shows maps of the /3lm parameters for all 
energy and time slices, along with the TRPES. This rep- 
resentation of the TRPADs is essentially one step further 
removed from the raw data, so again should be used in 
concert with plots showing the raw data before draw- 
ing firm conclusions, but also provides the most reduced 
and tractable form of the measurement. In this format, 
the 3 energy regions discussed in detail in the preceding 
can be readily compared. It is clear that the higher en- 
ergy region correlates to larger /? 2 ,o values, which change 
little over the main part of the signal, while the lower 
energy region shows marked oscillations; similar oscilla- 
tions extend to the lower energy slices, with peaks around 
t = 20, 70 fs. At later times, t > 60 fs, the trend is for 
reduced magnitude @2 y o, particular over the main signal 
in the 0.9 - 1.3 eV region. These times are outside of the 
cross-correlation of the laser pulses, so should be repre- 
sentative of only the parts of the excited state wavepacket 
which move orthogonal to the steep gradients on the po- 
tential energy surface responsible for the speed of the dy- 
namics and the rapid IP change observed (i.e. population 
which remains in a given region of configuration space 
for a longer time than the main part of the wavepacket). 
The /^o map shows much more oscillatory behaviour, 
which appears to show no obvious correlation with en- 
ergy or time, except for the trend towards larger positive 
values at t > 60 fs. For higher time resolution data it 
would be feasible to Fourier Transform this data to ex- 
tract the frequency content, but for the dataset shown 
here the limited number of delays results in only a crude 
frequency spectrum of little utility. 

In summary, the TRPADs presented here contain a 
plethora of information, and show very complex tempo- 
ral evolution, in contrast to the TRPES which provides 
little information on the temporal evolution of a given en- 
ergy slice over the probe pulse envelope, nor provides an 
observable as sensitive to the underlying dynamics. The 
challenge for the experimentalist is to determine whether 
the richness of the measured TRPADs can be interpreted 
in terms of the underlying dynamics without recourse to 
detailed theoretical treatments, that is to say without a 
priori knowledge of the underlying molecular dynamics 
or a full ab initio treatment of the dynamics and ioniza- 
tion. These points are discussed further below (section 
HIE). 



D. High-dimensionality mappings: correlated 
observables 



The advantages of measuring in coincidence have al- 
ready been discussed in terms of removal of background 
or other unwanted signals from the data hypercube. Nat- 
urally a further advantage is the ability to look for ion- 
electron correlations, and retrieve minor channels which 
would otherwise be inaccessible; for instance, the photo- 
electron spectra of the fragment channels, which would 
not be distinguishable in a non-coincidence measurement. 



Figure 10 shows the photoelectron spectra (time- 
integratedjfor the major ion channels as shown in fig- 
ure |5j The parent ion correlated spectra shows the main 
cross-correlation feature around 1 eV, and a rise towards 
lower energy, as already seen in the TRPES mapping (fig- 
ure [6]). As expected, the parent +1 channel has an identi- 
cal spectrum. However, the diffuse part of the parent ion 
signal has a very different spectrum, consistent with a 
different ionization process. Similarly, the fragmentation 
channels show different spectra; in all cases the spectrum 
is broad, and there is a drop or even disappearance of the 
1 eV peak seen in the parent ion channel. These spectra 
therefore provide additional information towards under- 
standing the excited state dynamics, with the spectra 
providing a fingerprint of different ionization channels. 
Here, the suppression of the cross-correlation peak in the 
fragment channels is consistent with the delayed opening 
of these channels, as observed in the time-resolved ion 
yields. 

Furthermore, the similarity between the channels sug- 
gests ionization occurs from a similar region of config- 
uration space in all cases. This is consistent with the 
similarity of the rise times and fall times observed in the 
time-resolved fragment yields, which intuitively suggests 
only a single dynamical pathway, shared by all the frag- 
ment channels, hinting at localization of the excited state 
wavepacket in configuration space. The similarity of the 
time-scales of the fragment channels, considered in light 
of the similarities in the correlated photoelectron spec- 
tra, therefore suggests that the fragmentation pathway 
of the ion is very sensitive to the form of the wavepacket 
at the time of ionization and, possibly, bifurcation on 
the ionic state (s) may lead to channels with apparently 
similar temporal response, but quite different fragmen- 
tation products. This is consistent with the picture of 
rapid wavepacket motion, with little dispersion, on the 
excited state and additionally indicates complex disso- 
ciation dynamics on the ionic surfaces. In terms of the 
experimentally accessed ionization pathways, the obser- 
vation of similar photoelectron spectra for the fragment 
channels is consistent with the possibility of sequential 
1+2' processes (ionization to Dq followed by absorption 



of a second probe photon) as discussed in section IIIB 
but suggests that direct 1+2' ionization to higher lying 
cationic states, and 1+1' ionization to Di, are unlikely 
channels. Direct ionization to excited ionic states would 
be expected to correlate generally with different spec- 
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FIG. 9. Energy-time maps of the total counts (TRPES) and extracted /3lm parameters. 




FIG. 10. Time-integrated photoelectron spectra correlated 
with the major mass channels (as shown in figure [5ja)). The 
parent ion correlated spectrum is scaled down by a factor of 
20 for plotting purposes. 



tra, although it is always possible that the spectra would 
not be significantly different - especially given the dif- 
fuse, unstructured nature of the photoelectron bands of 
butadiene - so this is not a rigorous conclusion. 

Various other correlated mappings are available from 
the data, and will be explored in a future publication. For 
example, the fragment correlated TRPES and TRPADs, 
although for weak channels the latter is very demanding 
statistically. For photodissociation studies energy corre- 
lation spectra - maps of electon vs. ion kinetic energy - 
are also useful 12 . 



E. Discussion - mapping dynamics with 
multi-dimensional measurements 



The data presented herein indicates that rapid and 
complex dynamics are present in butadiene, consistent 
with earlier experimental [38] [43] [44] and theoretical 
[SUES] works. The benefit of a multidimensional dataset 
is clear, with complementary information available from 
the coincident electron and ion data. The overall shape 
of the dynamics, that of a near- Gaussian wavepacket 
which moves along steep gradients on the potential en- 
ergy surface which are strongly correlated with the ver- 
tical IP, and shows little dispersion along other coordi- 
nates, emerges rapidly from the data. This picture fits 
both the chirped TRPES data and the delayed-onset of 
fragmentation channels which are assumed to be corre- 
lated with ionization to D\. The observation of TR- 
PADs which show very rich behaviour on rapid timescales 
(<20 fs) is striking, and indicates the sensitivity of the 
TRPADs to the dynamics under study. 

However, more specific mechanistic details, such as the 
nuclear motions involved or the mapping of these motions 
onto the TRPADs, do not seem to be forthcoming. This 
is in contrast to simpler cases, such as previous work on 
the NO dimer [16] or CS 2 [H [37], in which the lim- 
ited dimensionality of the problem enabled a more direct 
empirical approach, coupled with symmetry-based mod- 
elling to understand the measurements in more detail. 
Understanding the information conveyed by the TRPADs 
from larger systems therefore remains a significant chal- 
lenge. 

In the case of butadiene, detailed ab initio dynamics 
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calculations do exist [39], and further studies are ongo- 
ing [69]. It is therefore very tempting to compare the 
measurements to these calculations, despite the fact that 
these calculations do not, so far, explicitly include the 
ionization matrix elements. For instance, figures 5 and 6 
of ref. [39] show the behaviour and timescales associated 
with the bond alternation coordinate (a motion involv- 
ing bonds along the carbon backbone streching and com- 
pressing) and the out-of-plane twisting motions (twisting 
about the central CC bond and twisting of the terminal 
methyl groups). The bond alternation coordinate shows 
oscillations with a period of ~20 fs, and large amplitude 
motion with changes up to 0.4 A. The twisting coordi- 
nates show rapid changes over ~ 40 fs, followed by almost 
constant bond angles, with some fluctuations, for the re- 
maining 160 fs of the calculations. These motions are 
very suggestive of the oscillations observed in ft.oW f° r 
some of the low-energy slices, and the gradual rise over 
~60 fs observed in ^^(t) over all photoelectron ener- 
gies. However, from the data presented in ref. [39 it is 
not apparent how these motions map to the IP, which is a 
very significant factor in drawing these conclusions more 
firmly. As shown in ref. [29] such mappings may be very 
complex, and difficult to determine from anything but a 
full calculation including the ionization matrix elements. 

Another conclusion from the dynamics calculations, 
which is likely to be more robust, is the agreement be- 
tween the timescale of the population dynamics and the 
gross changes in both /?2,o(£) an d /?4,o(£) a t all energies 
for t > 60 fs. Figure 4 of ref. [39 shows the popu- 
lations of the bright and dark adiabatic states, defined 
there as Si and S2 respectively. The Si population dips 
almost immediately from its initial value at t = 0, with 
almost equal populations of Si and S2 at t « 25 fs, the Si 
population then grows again and peaks at around 50 fs, 
while the S2 population drops rapidly. Unsurprisingly, 
given the similarity of the timescales to those mentioned 
above, the non- adiabatic coupling between these states is 
primarily mediated by the out-of-plane twisting motions. 
In terms of the electronic character of the states, the 
early time dynamics is best considered in a diabatic pic- 
ture, with the bright state 1 B U character maintained as 
population switches from the Si to the S2 surface at early 
times [39], while over the longer timescale the electronic 
character becomes mixed by the out-of-plane twisting, 
and population is transferred non-adiabatically between 
62 and Si. 

Given that the adiabatic states are of different elec- 
tronic character one would expect a different PAD from 
each state; additionally, interferences between photo- 
electrons from the two states (different parts of a split 
wavepacket) may also play a significant role in the ob- 
servable if electrons of the same energy are created from 
the two states, and this could lead to more rapid mod- 
ulations in the form of the TRPADs (this is essentially 
the concept of the two-state model previously applied to 
CS2 [37 J, particularly in the case of symmetry-breaking 
leading to mixing of the adiabatic electronic state char- 



acters. With this in mind, it appears that the (adiabatic) 
population dynamics may have a more direct link to the 
observed TRPADs, with the changes at t > 60 fs corre- 
sponding to the region where most of the population is 
on a single adiabatic state, and the molecular geometry, 
hence electronic character, has stabilised after the ini- 
tial, rapid, out-of-plane twisting. Although this picture 
seems reasonable, without further modelling or calcula- 
tions these conclusions remain somewhat tentative. 

The challenge for the experimentalist therefore re- 
mains: to determine how the TRPADs, and other cor- 
related observables, can be interpreted without recourse 
to detailed theoretical treatments, hence without a pri- 
ori knowledge of the underlying molecular dynamics or 
a full ab initio treatment of the dynamics and ionization. 
While such computational approaches are very powerful, 
they are also time consuming, difficult, and only compu- 
tationally tractable for small molecules. Naturally the 
same applies to a purely experimental approach to the 
TRPES and TRMS data, but typically the lower dimen- 
sionality of this data lends itself to only broad interpre- 
tations, such as the presence of band switching [31] or 
quantum beats, which may fully describe the wavepacket 
dynamics in small systems [70 or arise from only the re- 
current parts of the wavepacket [18] [37] in more complex 
cases. Moreover, TRPES data is routinely interpreted 
via fitting of temporal profiles in what is essentially a 
principle component or single-value decomposition anal- 
ysis [TTJ [72] . While such analysis provides a quantitative 
breakdown of the data into spectral and temporal func- 
tions, it makes many assumptions regarding the under- 
lying dynamics (vis. that the observables equate directly 
to state populations in the standard interpretations, or 
at the very least behave phenomenologically in a statis- 
tical manner) and cannot be used to treat more complex 
dynamics, such as the kind observed here, without re- 
course to many time constants. The TRPADs clearly 
convey more information on the wavepacket dynamics, 
but the mapping from wavepacket to observable is inher- 
ently complicated. 

Aside from the TRPADs, the full imaging data pro- 
vides further correlated observables which may also be 
considered in light of the computational results. In 
particular, the time-resolved ion yields and fragment- 
correlated photoelectron spectra may provide a way of 
experimentally probing the three Si — >• So CIs discussed 
in ref. [39], which correspond to different molecular ge- 
ometries. It is likely that ionization from these different 
regions of configuration space (if energetically possible in 
the limited observation window) would lead to different 
photoelectron spectra and different fragmentation prod- 
ucts from D\. Even in the case where direct ionization 
from these regions is energetically forbidden, ionization 
of the parts of the wavepacket on Si heading towards 
these distinct regions of configuration space would also 
result in different dynamics on the ion surfaces, hence 
possibly lead to different fragmentation products. Fur- 
ther analysis of ion-electron correlations should present 
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deeper insight into the excited state dynamics than un- 
correlated ion or electron measurements alone. 



IV. CONCLUSIONS & FUTURE WORK 

The data presented here illustrates the richness of co- 
incidence imaging datasets, even in the case of very rapid 
and complex dynamics; it also highlights the difficulty of 
interpretation of such data. In this work, we have fo- 
cussed on the TRPADs obtained in coincidence with the 
parent ion channel, and made some tentative compar- 
isons with theory. In general the full 7D data provides 
the best hope of understanding the excited state dynam- 
ics of molecules but, unsurprisingly, significant effort is 
required to gain understanding directly from the data, 
or with the aid of theory. Despite the depth of material 
presented here there are still many avenues to explore 
in order to obtain more complete, unified experimental 
and theoretical approaches to problems in excited state 
dynamics. 

In future work we hope to continue bridging this gap 
between experiment and theory via a number of routes: 
by the more direct and explicit comparison of the exper- 
imental data with dynamics calculations, and the inclu- 
sion of observables in the calculations; by exploring fur- 
ther the possibilities of modelling for a qualitative/semi- 
quantitative treatment of the dynamics, expanding upon 
previous work which employed a simple 2-state model 
[37] . possibly to include a basic wavepacket treatment; 



through further analysis of the ion fragment distribu- 
tions, which have not yet been fully explored; by mak- 
ing further experimental measurements for aligned sys- 
tems and with different pump-probe polarisation geome- 
tries in order to obtain more detailed TRPADs (in both 
cases more L, M terms are allowed in the PADs), with 
harder VUV probe photons to access a larger photoelec- 
tron energy range, and with a stability enhanced set-up 
to allow for longer experimental runs; by further use of 
fragment correlated TRPES/TRPADs, so far not statis- 
tically feasible due to the dominance of the parent ion 
channel and consequent low total counts in these frag- 
ment channels. This multi-dimensional approach to a 
multi-dimensional challenge will hopefully prove fruitful 
for a deeper understanding of excited state dynamics in 
polyatomic molecules and, in particular, help to make 
the experimental measurement of such dynamics a more 
routine and insightful tool. 
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